function pOut = newparam4evolution(pBest,ncases,pOriginal,SDdivide)
% Based on best design, create more designs like it


% Initialize
pOut = [];
NameOfFields = fieldnames(pOriginal);
p = [];

for ii = 1:length(NameOfFields)
    temp = getfield(pOriginal,NameOfFields{ii});
    if numel(temp) == 1
        p = setfield(p,NameOfFields{ii},getfield(pBest,NameOfFields{ii}));
    elseif numel(temp) == 2
        p = setfield(p,NameOfFields{ii},[getfield(pBest,NameOfFields{ii}) temp(2)/SDdivide]);     
    else
        p = setfield(p,NameOfFields{ii},[getfield(pBest,NameOfFields{ii}) temp(2)/SDdivide temp(3)]);
    end
end

pOut = param4evolution(p,ncases-1);
pOut{ncases} = pBest;

